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1. INTRODUCTION 


The effects of atmospheric turbulence In both horizontal and 
near- horizontal flight, during the return of the Space Shuttle, are Impor- 
tant for determining design, control, and "pllot-ln-the-loop" effects. A 
non-recurslve model (based on von Karman spectra) for atmospheric turbu- 
lence along the flight path of the Shuttle Orblter has been developed which 
provides for simulation of Instantaneous vertical and horizontal gusts at 
the vehicle center-of-gravlty, and alro for simulation of Instantaneous 
gusts gradients. Based on this model the time series for both gusts and 
gust gradients have been generated and stored on a series of magnetic 
tapes which are entitled Shuttle Simulation Turbulence Tapes (SSTY). The 
time series are designed to represent atmospheric turbulence from ground 
level to an altitude of 120,000 meters. 

A description of the turbulence generation procedure Is provided 
in Section 2. The results of validating the simulated turbulence are 
described in Section 3. Conclusions and recommendations are presented In 
Section 4 with Section 5 containing references cited. Appendix A contains 
the tabulated one-dimensional von Karman spectra while Appendix B provides 
a discussion of the minimum frequency simulated. Appendices C and D present 
the results of spectral and statistical analyses of the SSTT. A more detailed 
description of the proper use of the tapes Is provided elsewhere [1]. 
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2. TURBULENCE GENERATION PROCEDURE 


The non-recurslve turbulence model used to generate the SSTT Is 
based on von Karman spectra with finite upper limits corresponding to the 
dimensions of the Space Shuttle, relative to the scale of turbulence In the 
atmosphere. Because the scale of turbulence Increases with altitude while 
the dimensions of the Space Shuttle are fixed, the finite upper limits of 
the von Karman spectra Increase with altitude. In order to take Into 
account the resulting spectral changes, the atmosphere, extending from 
ground level to 120,000 meters, was divided Into six altitude bands. The 
subsections which follow provide a description of the development and 
application of the turbulence generation procedures. 

2.1 BACKGROUND 

The current turbulence model represents the results of the 
development and evaluation of several different turbulence simulation 
techniques. Two of the earlier techniques, which were given serious 
consideration, warrant further discussion. 

Initial efforts involved refinement and evaluation of a turbulence 
model, TBMGD [2], which had been developed elsewhere [3]. This model was 
based on discretization of the Fourier Integral representation of turbulence 
and was designed for use with von Karman spectra. Several problems were 
encountered with TBMOD, both theoretical and practical. First, from a 
theoretical standpoint, the assumptions used in the development of the 
shear (gust gradient) simulation were difficult to justify. Second, from 
a practical standpoint, the output of TBMOD, representing turbulent gusts, 
when subjected to Fast Fourier Transform (FFT) spectral analysis, did not 
possess the proper von Karman spectral shape. Because of such problems, 
further development of the model was halted. 

The second simulation technique was based on digital filter 
theory coupled with a combined von Karman - Saffman spectral model (4) . 
Meromorphic Functions were used to approximate the various spectra. Based 
on z-transform theory, recursive difference equations were derived from 
such approximations. Such difference equations were then used to generate 


2-1 



the appropriate turbulence gusts and gust gradients. Unfortunately the 
recursive difference equations. In addition to being somewhat complex, 
proved numerically unstable and could not be used In their original 
recursive form [5]. 

A non-recursive version of the same difference equations was 
subsequently developed In an effort to overcome the stability problem [5]. 
These difference equations were also quite complex but resembled In form 
the output of digital filters, characterized by some Impulse response 
function, with a white noise Input. The present model, as described In 
subsections 2.2 through 2.4, was an outgrowth of this resemblance. 

2.2 SELECTION OF ATMOSPHERIC BANDS 

The standard deviations (oj, o 2 , o 3 ) and the Integral scale 
lengths (Lj, L 2 , L,) of atmospheric turbulence are functions of altitude, 
as shown in Table 2-1. Notice should be taken that the values for and 

presented In this table are consistent with those presented In JSC 7700 
{6]. Based on the variation of and presented In Table 2-1, the 
atmosphere was divided Into six altitude bands as presented In Table 2-2. 
Within each band, as also indicated in Table 2-2, characteristic Integral 
scales of turbulence were selected for use In calculating the finite 
upper limit of the turbulence spectral model discussed In subsection 2.3. 

2.3 DEVELOPMENT OF VON KARMAN SPECTRA WITH FINITE UPPER LIMITS 

As developed previously 15] the basic three-dimensional von Karman 
relation to be Integrated for the dimensionless gust spectra Is, 


^11 ^1 *^ 2 ’^ 3 ^ 


55 

36ai? 


(n 2 -^ 2 ) 

(1+Q 2 ) 17/6 


( 2 - 1 ) 


The corresponding von Karman relation for dimensionless gust gradient 
spectra is 




55 

*77 (i*7) T77r 


( 2 - 2 ) 


These tnree-dlmenslonal spectral relations must be Integrated over certain 
ranges of values of t? 2 and 0 3 to obtain one-dimensional spectral models 
and ♦n/jj( n 1 )- 
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TABLE 2-1. VARIATION OF STANDARD DEVIATION 
AND LENGTH SCALE WITH ALTITUDE* 



Double entries for a tabulated altitude Indicate a step 
change In standard deviation or Integral scale at that altitude. 
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TABLE 2-1. VARIATION OF STANDARD DEVIATION 

AND LENGTH SCALE WITH ALTITUDE (Continued) 


ALTITUDE 

(m) 

STANDARO DEVIATION 
OF TURBULENCE 

INTEGRAL SCALES 
OF TURBULENCE 

a l (m/sec) 

a 2(m/sec) 

°3(m/sec) 

Lj(ni) 

L 2 (m) ! 

i 3 («i 

27000 

7.00 

7.00 

4.22 

20000 

20000 

1230 

30000 

8.23 

8.23 

4.66 

23533 

23533 

1443 

40000 

12.82 

12.82 

6.09 

36693 

36693 

2231 

50000 

18.08 

18.08 

7.51 

51786 

51786 

3128 

60000 

23.94 

23.94 

8.90 

68623 

68623 

4124 

70000 

30.36 

30.36 

10.28 

87063 

87063 

5208 

80000 

37.29 

37.29 

11.65 

106998 

106998 

6376 

90000 

44.70 

44.70 

13.01 

128338 

128338 

7622 j 

100000 

52.58 

52.58 

14.35 

151010 

151010 

8941 

110000 

60.89 

60.89 

15.69 

17495C 

174950 

10330 

120000 

69.62 

69.62 

17.02 

200000 

200000 

lieoo 


2.3.1 Upper Limits of Integration 

The upper limits of Integration for j = 2 and 3 are calculated 
according to the relation [7] 

n 1jmax * aL 1 /l j lJ * 2 - 3 > U ' 3) 

where 

a = 1.339 

L. 3 integral scale of turbulence associated with the 
1 ( n l) s P* ctrum 

l. * characteristic length of Space Shuttle In the 
J jth direction 

Values of L^ for the six bands are given In Table 2-2 while the character- 
istic lengths, * j . are presented In Table 2-3. 
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TABLE 2-1. SUMMARY OF TURBULENCE PARAMETERS 
IN DISCRETE ALTITUDE BANDS 
























































TABLE 2-3. CHARACTERISTIC DIMENSIONS 
OF THE SPACE SHUTTLE [8] 


Characteristic Magnitude 

Length (ft) (m) Explanation 


*1 

39.56 

12.06 

mean aerodynamic chord 

£ 2 

39.05 

11.9 

1/2 wingspan 

£ 3 

10.95 

3.34 

1/2 fuselage thickness 


In the case of special consideration must be given to the 

dimensional frequencies corresponding to the dimensionless limits. The dimen 
slonal frequency limit satisfies the relation 

^lmax 3 n 11*x V '< 2w "-1> < 2 ' 4 > 

where 

V = vehicle velocity 

The maximum dimensional frequency which the Space Shuttle simulators are 
capable of handling is 4 hertz. Thus any higher frequencies should be ex- 
cluded from tho simulation. For this reason the dimensionless frequency 
limit must satisfy the relation 

n Unw>x ' min < al V t l- 2ML 1 f lmax /V> < 2 ' 5 > 

where 

f lmx * 4hertz 

Values of ^ lmax based on Eq (2-5) are Included In Table 2-2. 

2.3.2 One-Dimensional Spectra 

There are six spectra of primary Interest for turbulence simu- 
lation, as indicated in Table 2-4. Based on second-order numerical Inte- 
gration, the six corresponding three-dimensional gust and gust gradient 
spectral relations, as given by Eqs (2-1) and (2-2), were Integrated over 
Oj and ^ (with the appropriate upper limits). The resulting one-dimensional 
spectra for all altitude bands are presented In Appendix A. These spectra 
were used In establishing the Impulse response functions associated with 
digital filter simulation processes described In subsection 2.4. 
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TABLE 2-4. TYPES OF SIMULATED TURBULENCE 


Type 

Corresponding 

Spectrum 

Comments 

Jt 

*n 

longitudinal gust 

u 2 

*22 

transverse gust 

u 3 

*33 

vertical gust 

dU2/dXj 

*22/11 

yaw 

9u-j/3Xj 

*33/11 

pitch 

3u 3 /3x2 

*33/22 

roll 


2.3.3 Dimensionless Energy Content 

The total dimensionless energy content of each one-dimensional 
spectra In each altitude band was established by integrating the corres- 
ponding spectra over the appropriate finite limits, indicated in Table 2-2. 
The resulting energy content is presented In Table 2-5. As might be expected 

TABLE 2-5. DIMENSIONLESS ENERGY CONTENT 
FOR GUSTS AND GUST GRADIENTS 


ALTITUDE 

SPECTRUM 

BAND 

*11 

*22 

‘33 

*22/11 

*33/11 

*33/22 

l 

.6225 

.5010 

.2752 

.5877 

.1525 

.1557 

2 

.8595 

.8560 

.8383 

13.147 

12.171 

12.308 

3 

.8956 

.8952 

.8809 

24.767 

22.643 

22.890 

4 

.9298 

.9296 

.9197 

54.123 

49.527 

50.060 

5 

.9977 

.9953 

.9251 

1740. 

41.71 

95.62 

6 

1.000 

.9973 

.9363 

2309. 

52.08 

391.6 
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the total dimensionless energy content of each of the turbulent gust series 

* 

Is less than unity. The dimensionless energy content for each gust gradient, 
however. Is not limited In such a manner and range as high has 391.6. For 
both gusts and gust gradients the total energy content Increases with altitude 
because of similar Increases In the limits of Integration. 

2.4 DIGITAL FILTER SIMULATION 

As suggested In subsection 2.1, simulated turbulence, Y(t), can be 
Interpreted as the response or output of a control system [9] with double- 
sided response functions, h(t), subject to an Input consisting of Gaussian 
white noise I (t ) . This response can be represented by the convolution 
Integral 

Y(t) - l a h(x) I(t-x) dx (2-6) 

Based on filter theory the double-sided spectrum, ♦ D y(fi), of the simulated 
turbulence satisfies the relation 

♦dy^) = (2-7) 

where H(flj) = F{h(t)] 

^DI (^1 ) = double-sided power spectrum for white noise 

Generally the standard deviation of any white noise signal has a value of 

unity. Furthermore in most practical situations the white noise is defined 

** 

to occur over some Interval extending from -Q, to +n, . For this case , 

I max lmax 

* Di(ni) ■ ^ 

* N ( 2 - 8 ) 

f “ 

Actually the term "energy" is not precise when dealing with gust gradients. 

**The subscript, 1, normally applied to the variables n, and T, has been 
suppressed In Eqs (2-8) through (2-13) for simplicity. 
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where 


By substitution, 
♦qy^i) 


time Interval associated with generation 
process (-w/O^) 


H (ni)H*(«i) ^ 


(2-9) 


If H(flj) Is limited to real values, 

♦o»< a l> - H 2 ( n i>s 

Rearrangement of Eq (2-10) yields 

h(«i) ■ VrV^J 


( 2 - 10 ) 


( 2 - 11 ) 


Then based on the definition of the Inverse Fourier transform, the double- 
sided impulse response function h(t) can be expressed as 


,-l, 


h(t) - F IH(flj)] 


00 

./V^T cost^t) dflj 
•00 

00 

* 2 /^T^OV (n l) «0*(«1*) dn l 


cosfft^t) dflj 


♦yfflj) COS(fljt) dflj 


(2-12) 


where 




single-sided spectrum of Y(t) 


2-9 


The single-sided spectra tabulated In Appendix A correspond to ^y(n^) 

The discrete version of the convolution Integral given In Eq (2-6) yields 


Y(k) « £ h(j) I(k-j) T 
j-N 


(2-13) 


where 


Y(k) ■ discrete sampled turbulence outpu. 

h(j) » discrete double-sided Impulse response 
function h(jT) 

I(k) ■ discrete sampled white noise Input 


Eq (2-13) represents the basic, non-recursive relation for the generation of 

* 

simulated turbulence . The impulse response functions h(t) were evaluated 
by means of second-order numerical integration of Eq (2-12) using the six 
spectra from Appendix A. In carrying out this evaluation some maximum value of 
t must be established, corresponding to the value of N for Eq (2-13), and re- 
placing the infinite limit of Eq (2-6). As discussed In Appendix B this maximum 
time limit determines the minimum frequency, n^ lmax , for which the corresponding 
spectrum is accurately simulated. 

The values of dimensionless time Increment, T. , used for the six 
altitude bands are included in Table 2-2 and are based on the values of n^ max 
shown in the same table. Thus the Nyquist generation frequencies for the 
simulated turbulence correspond to the upper frequency limits for 0^ as com- 
puted by Eq (2-5) for each altitude band. 


In certain references 19,11] to correct for the "effect of digitizing" the 
series represented by Eq (2-13) has been divided by A. This process can be 
seen to be dimensionally Incorrect and actually results from the use of a white 
noise spectrum with unit strength Instead of a strength of T/2 tt. 
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2.5 


EFFECTS OF DIGITIZATION 


The effects of digitization In turbulence simulation have been 
considered by a number of Investigators [9-12]. As a result of these studies 
two basic digitization effects have been generally Identified. 

The first effect results from the assumption of a white noise 
specrrum with unit strength Instead of unit power, noted In subsection 2.4. 

To correct for such an "effect" the proposed procedure Is to divide the 
series approximation of the convolution Integral by A. This "effect" 
disappears when the white noise spectrum has unit power. 

The second effect involves the tapering of the spectrum of simulated 
white noise, <Pj , in the vicinity of the Nyqulst generation frequency, 

Some investigators [10,12] have considered It necessary, because of the taper- 
ing effect, to generate the simulated turbulence time series at a rate from 
four to ten times the rate at which the series will be sampled. 

The second effect arises from the discrete processes associated 
with both the generation and sampling of the simulated turbulence. In the 
case of discrete white noise with unit variance, the time series involved is 
basically a train of step functions as shown in Figure 2-1. The autocorre- 
lation function of the train of step functions depicted in Figure 2-1 can 



2-11 


be shown to be 


V T) 



(M ± t g) 

(|t| >t g ) 


The corresponding double-sided power spectrum by definition is 

♦' DI (fl) s F(R di (t)1 

T q sin 2 (nT Q /2) 

27 («T g /2) 2 

j s1n 2 (Oir/2n NG ) 

^ (fttr/2n NG ) 2 

The single-sided version of this power spectrum is shown in Figure 2-2 



Figure 2-2. White Noise Spectra 


(2-14) 


(2-15) 
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Theoretical white noise, by definition. Is characterized by a 
uniform power spectral distribution. In order to avoid Infinite power, such 
a spectral distribution Is normally restricted to the frequency band 
( -n NG - n - n NG^ * or a double " s * ded spectrum. For a signal with unit power 
the spectral density function for such white noise Is 


( - n NGl n ^ fl NG ) 

l o |n|) 


(2-16) 


Such a theoretical distribution In single-sided form Is also shown in Figure 
2 - 2 . 

It Is important to note that the two power spectra shown in Figure 
2-2 are both normalized and thus 


i 

f 


♦ DI (n) d « 


/ 




(2-17) 


- 1 

The theoretical spectrum Is basically a rectangular pulse function while the 
discrete spectrum Is characterized by tapering. The difference between these 
two spectra is generally considered the basis for the second digitization 
effect. 

The preceding descriptions of the two spectra <t> j (n) and 4>gj(fl) are 
based purely on mathematical theory. To observe such spectra in reality the 
corresponding time series would have to be sampled with an infinitesimal 
sampling Interval. Actually, finite sampling intervals, T $ , must be used but 
this finite (or discrete) sampling process results In aliasing. The aliased 
spectrum, ^(n), based on the finite sampling process, Is related to the 
original spectrum according to the relation [13] 
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(2-18) 


where 


♦♦(n) • E ♦«H2kn NS ) 
k *-«• 

fl^S ■ Nyquist sampling frequency (-ir/T^) 


For the case of the discrete white noise 

■ £ < ' 0I ( n * 2k! W 

k— 



sin 2 [(n+2kn NS )Tr/2n NG l 

[(n^2kn NS )TT/2P NG j 2 


(2-19) 


Numerical evaluation or this series has been carried out for ■ 100 with 
various ratios cf Including .5, 1, 2, and 4. The resulting aliased 
spectra are presented In Figure 2-3. It Is important to note that the figure 
Indicates that 


/ * 
^NG 


♦•V“> 


l-^G - a - "MG* \ 

/ ^HG * "NS* 

(°f(G < l n D 1 


( 2 - 20 ) 


In this case, by comparison with $ DI (0), 

*’* 01 tn) * ("ng'^ns’ (2 ' 2 

Thus for white noise the aliasing due to discrete sampling exact. y offsets 
the tapering due to discrete generation when the sampling frequency equals 
the generation frequency. Based on this fundamental point, it Is clear that 
In the simulation of white noise no tapering of the spectrum occurs as long 
as the sampling rate equals the generation rate. Under most conditions this 
equality Is automatically satisfied. 
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The process of convolving the white noise with the appropriate 
Impulse response function Is also carried out In a discrete manner. The 
process Involves selecting (or sampling) values of the white noise signal 
and the Impulse response function at equal Intervals In time and then approxi- 
mating the convolution Integral by a sui.mat1on of products. It Is Important 
to note that the discrete sampling of both the white noise and the Impulse 
response function normally occurs at the sam: rate as the generation rate 
for the white noise. Thus the resulting spectra for the sampled discrete 
white noise, as previously shown, will be uniform. According to the con- 
volution theorem, the spectrum of the output signal equals the spectrum of 
the Input white noise multiplied by the product of the Fourier transform of 
the impulse response function and Its complex conjugate. Thus, as previously 
noted in subsection 2.4, for a continuous signal, 

* DY (n) = ♦ DI (n)H(n)H*(n) ( 2 - 22 ) 

The Fourier transform H(£l) for the continuous impulse function, h(t). Is 
H(fl) = F[h(t)] 

= VrV fl > (2 ' 23) 

The corresponding output spectrum for a discrete signal would be 

*'dy («) a ♦ , DI (n)H'(n)H'*(n) (2-24) 

The Fourier transform H'(n) of the discrete impulse response function, h'(t), 
Is 

H'(n) = F[h* ( t ) ] (2-25) 

Based on the preceding development, for cases In which the sampling 
frequency equals the generation frequency, any difference between the discrete 
turbulence spectrum and the continuous spectrum apparently originates because 
of some difference between H(fl) and H’(n). 
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WHITE NOISE ENERGY DENSITY 



Figure 2-3. Effects of Aliasing on White Nose Spectrum 
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3. SIMULATED TURBULENCE TAPES 


The turbulence generation procedure described In Section 2 has been 
used to generate six dimensionless simulated turbulence time series which are 
stored on magnetic tapes as summarized In Table 3-1. The appropriate procedures 
for using the tapes are described elsewhere [1]. Subsection 3.1 provides a 
description of the results of validating the tapes while subsection 3.2 presents 
an explanation of the process for converting from dimensionless to dimensional 
values. 


Tape 

TABLE 3-1. INDEX OF SHUTTLE SIMULATED 
TURBULENCE TAPES (SSTT) 

Time Series 

Comments 

SSTT-1 

u^ - gust 

longitudinal gust 

SSTT-2 

u 2 - gust 

transverse gust 

SSTT-3 

u^ - gust 

vertical gust 

SSTT-4 

3u2/3Xj - gust gradient 

yaw 

SSTT-5 

BUj/dXj - gust gradient 

pitch 

SSTT-6 

3u 3 /3x ? - gust gradient 

roll 

3.1 

VALDIATION OF SIMULATED TURBULENCE 



A spectral analysis of each of the dimensionless time series has been 
carried out by means of a Fast Fourier Transform FFT4 [141 . The results, which 
are presented In Appendix C, demonstrate that the simulated turbulence possesses 
the proper von Karman spectral characteristics. 

All of the dimensionless time series have also been analyzed statis- 
tically to determine the gust and gust gradient probability density functions. 

As shown In Appendix D the results of these analyses Indicate that both the 
simulated gusts and gust gradients are normally distributed, with near-zero 
means and standard deviations consistent with the energy content presented 
in Table 2-5. 
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3.2 ' CONVERSION TO DIMENSIONAL VALUES 

The dimensionless time series on each tape must be converted to 
dimensional form before actual use In a simulation exercise. The conversion 
process generally Involves multiplication anl/or division by the appropriate 
turbulence parameters. For dimensionless gusts, u^, the corresponding standard 
deviation, o^, should be used. Thus 

u* « o i u 1 (3-1) 

where 


* 

Uj ■ dimensional gust 


For dimensionless 
Thus 


3Uj 


gust gradient, , the parameters and Lj are 

J 


used. 


♦ 

3u . 

o, 3u, 


X — — 

3Xj 

L 1 sx j 

where 


* * 


3Uj/3Xj 

■ dimensional gust gradient 


(3-2) 


In the case of dimensionless time It Is necessary to develop the 
procedures for converting both from dimensionless to dimensional form anJ also 
to dimensionless from dimensional. In proceeding from dimensionless to dimen- 
sional time the dimensionless time step, T^, represents the basic unit to be 
converted. The conversion Involves the vehicle velocity, V, and the turbulence 
scale, Lj. Thus 

At* - aL^/V (3-3) 


where 



dimensional time step 
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It Is Important to nqte that because both Lj and V vary with altitude, 

the resulting dimensional time step At^ Is not a constant. To obtain dimensional 
* 

time, t^, a summation process Is Involved as follows: 


where 



N 

= E At 
n=l 


* 

In 


N 

aT ( E 

Vl 


L 1„' V n 


(3-4) 



w 

V(Z„) 

altitude at nth step 


In converting to dimensionless from dimensional time the basic 

★ 

unit, the dimensional time step, fit , will normally be a constant. The 
corresponding dinrrsionless time interval, T- m , will be 


T. 

im 


V fit 

m 


* 



(3-5) 


The total dimensionless time, t^, will be 


t. 


iW 


M 


H V_it 

E n 


m=l aL im 

* M 

^£V L < 


m 


(3-6) 
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The dimensionless time, t^, corresponds to some M' dimensionless time Inter- 
vals, Tj , plus some fractional interval, T', as follows: 

*m ■ M ' T i * T ' (3-7) 

where 

0 < T' < T i 


ihus the number of dimensionless time intervals, M', can be computed 
as 

M 1 = Int(t^ M /T 1 ) 

= Int{ aT7 m ?-| V m /L 1m^ 

i iii s 1 

where 

Int ( ) = integer value of ( ) 

The fractional interval, T' , can be computed by the relation 
r = t iM -M'T. 


(3-8) 


(3-9) 


The interpolation process will involve interpolating between t.^, and t^,^ 
at the point t iM as shown in Figure 3-1, 




t iM'-l 1 1M‘ f t iM‘>l 


Oimensionless 

Time 


'1H 


Figure 3-1. Relationship Between t^, t^,, and t^.+j 


In the conversion to or from dimensional values three parameters 
are required: standard deviation, integral scale of turbulence, and vehicle 
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speed. The variation of the turbulence standard deviation, a^, with 
altitude was presented In Table 2-1. The same table contains the turbulence 
scale, Lj, as a function of altitude. The vehicle speed, V, Is a function 
of altitude but also may vary from one trajectory to another. Table 3 -2 
provides representative values of V as a function of altitude. 

TABLE 3-2. VARIATION OF SHUTTLE SPEED 
WITH ALTITUDE [12] 


ALTITUDE 

(m) 


100 


300 

m 

500 

158 

2000 

170 

4000 

188 

6000 

200 

8000 

240 

10000 

300 

20000 

500 

40000 

1928 

60000 

4695 

80000 

7468 

100000 

7521 

120000 

7510 





4. CONCLUSIONS AND RECOMMENDATIONS 


By means of a non-recurslve discrete generation process, based on 
a von Karman spectral model with finite upper limits, dimensionless simulated 
turbulence time series have been developed and stored on six magnetic tapes. 
Longitudinal, transverse, and vertical gusts are simulated as well as the 
gust gradients associated with yaw, pitch, and roll. For each gust or gust 
gradient six separate time series (corresponding to the six altitude bands 
extending from ground level to 120,000 meters) have been stored on each tape. 

The results of spectral analyses of each tape reveals that the 
simulated turbulence possesses the appropriate von Karman spectral character- 
istics. Statistical analyses of the tapes Indicate that both the simulated 
gust and gust gradients are normally distributed with near-zero means. Further- 
more the standard deviation of each series is constant with the theoretical 
energy content. 

The Shuttle Simulated Turbulence Tapes (SSTT) are now ready for 
actual use for simulating turbulence at altitudes up to 120,000 meters. 
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APPENDIX A 


DIMENSIONLESS VON KARMAN SPECTRA 
WITH FINITE UPPER LIMITS 


For each altitude band, the three-dimensional spectral model 
for gusts, as given by Eq (2-1), and the three-dimensional model for gust 
gradients, as given by Eq (2-2), have been Integrated with respect to ^ 
and Oj over the finite limits calculated according to Eq (2-3). The six 
resulting one-dimensional spectra are presented in Tables A-l through A-6, 
corresponding to Altitude Bands #1 through 6 respectively. These spectra 
were used in the numerical evaluation of the Impulse response functions 
described in subsection 2.4. 
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APPENDIX B 


ESTABLISHMENT OF LOWER FREQUENCY LIMITS 


The maximum time limit, t max , for which the Impulse response 
function Is computed, determines the minimum frequency, for which 

the corresponding spectrum Is accurately simulated according to the relation 


fi 11 min a "^max 


(B-l) 


The simulated turbulence may contain lower frequencies, depending on the total 
length of the time series, but the shape of the actual spectrum for any fre- 
quencies less than ^ 1(n ^ n will not In general match the theoretical spectral 
shape. To verify this point for the $ 33/11 spectrum two impulse response 
functions were generated, the first extending out to a t max of 42.10573 and 
the second extending out to a t max of 85.895689. Two separate turbulence 
time series were then generated, one for each impulse function. The results 
of spectral analysis of these two time series Is presented in Figures B-l 
and B-2. As Indicated In each figure, the observed spectrum tends to drift 
away from the theoretical spectrum for frequencies below ir/ , t inax . The impulse 
response function with the larger t produces a time series whose spectral 
shape matches the theoretical spectrum to the lower frequency. 

Based on various characteristics of the Space Shuttle simulators, 
a minimum frequency, fi min » of .04 hertz was established for generating 
simulated turbulence. The dimensional frequency, fj min , is related to the 
dimensionless frequency, i m ^ n , according to the relation 

°1 lmln ' 2ML 1 < B ' 2 i 

By substitution and rearrangement, 

■ v '< 2sL i <W < B - 3 > 


To satisfy the requirement for a minimum simulation frequency of .04 hertz, 
values of each impulse response function were computed for 100 dimensionless 
time intervals. As shown in Table B-l, by using a constant number of time 


B-l 






Intervals, the actual minimum simulation frequency varies from .0174 hertz 
In Band 1 to .04 hertz In Band 6. Thus, at the lower altitudes the turbulence 
simulation will actually be valid for frequencies somewhat lower than the 
required minimum as given by Eq (8-2). 


TABLE B-l. DIMENSIONAL AND DIMENSIONLESS MINIMUM FREQUENCY LIMITS 


MINIMUM SIMULATION 
FREQUENCY LIMITS 


1 

SPECTRA 

1 

*11 

2 

*22’ *22/11 

3 

*33* *33/11 
’*33/22 


DIMENSIONAL 


.0174 

.0212 

.0211 

.0396 

.04 

.04 


174 

206 

211 

396 


.04 

.04 


.0174 

.0208 

.0211 

.0396 

.04 

.04 


DIMENSIONLESS 




11ml n 


.03075 

.2110 

.3331 

.5918 

7.365 

8.949 


.01865 

.2132 

.3331 

.5918 

.4630 

.5280 
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J T APPENDIX C 

I SPECTRAL ANALYSIS OF SIMULATED TURBULENCE 

By means of a Fast Fourier Transform (14] spectral analyses of all 
1 * * 

simulated turbulence have been performed . The results are presented In 

dimensionless form In Figures C-l through C-36. Table C-l provides a summary 

of these figures. Also Included In each figure Is the theoretical von Karman 

spectra. The agreemcw between the theoretical spectra and the computed spectra 

Is quite satisfactory. 


TABLE C-l. MATRIX OF SPECTRAL ANALYSIS FIGURES 


SERIES 

TYPE 

ALTITUDE BAND 

1 

2 

3 

4 

5 

6 

U 1 

C-l 

C-2 

C-3 

C-4 

C-5 

C-6 

u 2 

C-7 

C-8 

C-9 

C-10 

C-ll 

C-12 

U 3 

C- 13 

C- 14 

C- 15 

C- 16 

C- 17 

C-18 

Su^/Bx j 

C- 19 

C-20 

C-21 

C-22 

C-23 

C-24 

3u 3 /3Xj 

C-25 

C-26 

C-27 

C-28 

C-29 

C-30 

3u 3 /3X2 

C-31 

C-32 

C-33 

C-34 

C-35 

C-36 



The spectral analysis Involved the first 4096 terms of each time series except 
for bands 5 and 6 for the u^ and U£ gusts. For these cases 8192 terms were used. 




C-l 




Figure C-l. Uj - Gust Spectrum, Altitude Band #1 Figure C-2. Uj - Gust Spectrum, Altitude Band #2 


THEORETICAL THEORETICAL 

VON KARMAN VON KARMAN 

\ OBSERVED \ OBSERVED 




























Spectrum, Altitude Band #3 Figure C-16. u, - Gust Spectrum, Altitude Band #4 
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Gust Gradient Spectrum, Figure C-22. 3 u 2 /3Xj - Gust Gradient Spec*i 

and *3 Altitude Band #4 












OBSERVED 



3u 3 /axj - Gust Gradient Spectrum, Figure C-28. au-j/axj - Gust Gradient Spectrian 

Altitude Band #3 Altitude Band #4 
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APPENDIX D 


STATISTICAL ANALYSIS OF SIMULATED TURBULENCE 

By means of standard statistical analysis procedures each of the 
SSTT has been analyzed to determine its mean value, standard deviation, and 
probability density distribution . The resulting mean values are presented 
in Table D-l while Table D-2 contains the resulting standard deviations. 

As expected all mean values were near zero. The standard deviations represent 
the square root of the energy content. The rutio of the theoretical energy 
content (from Table 2-5) to the square of the corresponding standard deviation 
(from Table D-2) is presented in Table D-3. The agreement appears quite satis- 
factory. 

The gust and gust gradient probability density distributions are 
presented in Figures D-l through D-36 in accordance with Table D-4. In each 
figure the corresponding theoretical normal distribution is also presented. 

The results indicate that both the gust and gust gradient time series are very 
dost to normal distributions. 

TABLE D-l. MEAN VALUE OF GUST AND GUST GRADIENTS 


SERIES 

ALTITUDE BAND 

TYPE 

1 

2 

3 

4 

5 

6 

U 1 

-.019295 

-.042050 

-.051852 

-.083142 

-.0455 

-.0464 

u 2 

-.010671 

-.029576 

-.0371 

-.049431 

-.0428 

o 

ft 

u 3 

-.006U06 

-.029652 

-.037043 

-.049370 

-.043788 

-.046637 

i^/axj 

-.000002 

-.001572 

-.002794 

-.005628 

-.172152 

-.206385 


-.000001 

-.001591 

-.002798 

-.005649 

-.004293 

-.004893 

au^/axj 

-.005760 

-.073072 

-.103823 

-.160303 

-.171448 

-.288178 


The statistical analysis involved the first 4096 terms of each time series except 
fcr bands 5 and 6 for the and U£ gusts. For these cases 8192 terms were used. 
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TABLE 0-2. STANDARD DEVIATION OF GUST AND GUST GRADIENTS 


SERIES 

ALTITUDE BAND 

TYPE 

1 

2 

3 

4 

5 

6 

U 1 

.788959 

.927098 

.946351 

.964271 

.99888 

.99996 

u 2 

.707845 

.925201 

.94619 

.964152 

.99764 

.99863 ; 

u 3 

.524571 

.915606 

.938552 

.958985 

.961845 

.967651 

3uo/axj 

.766627 

3.625937 

4.976618 

7.356808 

41.717292 

48.052734 


.390512 

3.488677 

4.758426 

7.037516 

6.458092 

7.216648 


.394539 

3.508280 

4.784378 

7.075349 

9.778736 

19.790119 




TABLE D-3. RATIO OF THE THEORETICAL ENERGY CONTENT 

TO THE SQUARE OF THE OBSERVED STANDARD DEVIATION^ 


I 

| i 

\ l 

| 



■11 

ALTITUDE BAND 

1 

2 

3 

4 

5 

6 

U 1 

1.0001 

1.0000 

1.0000 

1.0000 

.9999 

1.0001 

u 2 

.9999 

1.0000 

.9999 

1 . 0000 

1.0000 

1.0000 

u 3 

1.0001 

1.0000 

1.0000 

1.0001 

1.0000 

.9999 i 

ai^/axj 

1.0000 

1.0000 

1.0000 

1.0000 

.9998 

1.0000 

3u 3 /3x 2 

1.0000 

1.0000 

1.0000 

1.0000 

1.0001 

1.0000 

auj/axj 

1.0002 

1.0000 

1.0000 

1.0000 

1.0000 

.9999 
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Figure D-16. u~ - Gust Probability Density Distribution, Altitude Band #4 
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